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We consider the energy spectrum and the spin-parity structure of the eigenstates for a quantum 
dot made of a strong topological insulator. Using the effective low-energy theory in a finite-length 
cylinder geometry, numerical calculations show that even at the lowest energy scales, the spin 
direction in a topologically protected surface mode is not locked to the surface. We find "zero- 
momentum" modes, and subgap states localized near the "caps" of the dot. Both the energy 
spectrum and the spin texture of the eigenstates are basically reproduced from an analytical surface 
Dirac fermion description. Our results are compared to microscopic calculations using a tight- 
binding model for a strong topological insulator in a finite-length nanowire geometry. 

PACS numbers: 73.50.-h, 72.80.Vp, 73.23.-b 



I. INTRODUCTION 

The theoretical prediction and subsequent experimen- 
tal verification of the conducting surface state of a strong 
topological insulator (TI) continues to generate a lot of 
excitement in physics; for reviews, see Refs.llHSl In a TI, 
strong spin-orbit couplings and band inversion conspire 
to produce a time-reversal invariant topological state dif- 
ferent from a conventional band insulator. Using Bi2Se3 
as a weakly correlated reference TI material with rather 
large bulk gap A;, ss 0.3 meV, surface probe experiments 
have provided strong evidence for the topologically pro- 
tected gapless surface state.!^ The measured spin texture 
of the surface state is consistent with predictions ob- 
tained for two-dimensional (2D) massless Dirac fermions. 
Under this "relativistic" description, the spin direction 
is locked to the surface, and the surface state is stable 
against the effects of weak disorder and weak interac- 
tions (topological protection) .EH^l is thus useful to first 
study the simplest case of a noninteracting disorder-free 
model, which is the case investigated below. 

Because of the residual bulk conductivity of the 
presently available (nominally insulating) TI samples, it 
has been difficult to experimentally extract the surface 
contribution to the electrical conductivity. One attempt 
to improve the situation is to consider mesoscopic sam- 
ples, where the surface-to-volume ratio is more advanta- 
geous. In particular, thin- film geometric and quasi- 
ID nanowires ( "ribbons" ^pESl have been studied exper- 
imentally. Signatures for Aharonov-Bohm interference 
effects associated with the topological surface state in 
Bi2Se3 nanowires were reported,^ cf. also related experi- 
ments for Sb2Te3 nanowires.l^ For infinitely long and cir- 
culary symmetric topological nanowires, band structure 
calculations predict a multi-channel waveguide where 
all surface modes are gapped because of spin-surface 
lo eking .EiHIll With k denoting the conserved momentum 



along the wire axis (taken along the z direction), and j 
the half-integer total angular momentum, the dispersion 
relation of these modes is 

E,,±{k) = ±^{v,kf + {]V2/R)\ (1.1) 

where h = 1 throughout and ± for conduction and va- 
lence band, respectively. The Fermi velocities vi and V2 
differ because the bulk dispersion relation is anisotropic, 
see below, and the nanowire radius is R. Note that there 
is a minimal gap = V2/R for the surface modes since 
j is half-integer. For reasonable values of R, we have 
A, < A,,. 

Experiments probing quantum dot physics in finite- 
length TI nanowires are expected to yield new insights 
into the exciting physics of TIs, in close analogy to semi- 
conductor nanowires and carbon nanot ubes where such 
experiments have been highly successful .I^^^H To prepare 
the ground, we here consider the band structure for TI 
quantum dots with a finite-length nanowire geometry. 
For a spherical TI dot, the band structure was worked 
out hefore^^^^ However, we draw attention to several 
features that only arise when the surface contains sharp 
edges, i.e., non-differentiable parts, as is the case for the 
cylinder. We employ three different and independent ap- 
proaches to understand TI quantum dot energy levels 
and their spin texture: (1) For a cylindrical TI nanowire 
of length L and radius R closed by flat caps, we have 
performed detailed numerical calculations for the energy 
spectrum and the spin texture of the eigenstates based 
on the effective low-energy theory of Zhang et Q/, po | 2i | 
material parameters were chosen for Bi2Se3 as quoted 
in Ref. 21, (2) An analytical approach starting from a 
surface Dirac fermion description has been developed for 
the same geometry. Most of our numerical results can 
thereby be quantitatively reproduced within an analyti- 
cal theory. (3) We have also studied a microscopic tight- 
binding model for a strong TI in the finite-size nanowire 
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Ao (eVA) 


Bo (eVA) 


Co (eV) 


Ci (eVA") 


Ca (eVA") 


Mo (eV) 


Ml (eVA^) 


Ma (eVA^) 


3.33 


2.26 


-0.0083 


5.74 


30.4 


-0.28 


6.86 


44.5 



TABLE I: Parameter values in Eq. \2.1\ appropriate for 
BiaSes (taken from Ref. l2T|l . 



geometry and find qualitatively similar results. 

The main conclusions reached from these three ap- 
proaches are as follows: First, longitudinal momentum 
quantization implies a discrete sequence of energy lev- 
els which can (roughly) be approximated by letting k — >■ 
kn = mr/L {n E N) in Eq. (1.1). Remarkably, in ad- 
dition we find unconventional "zero-momentum" states 
(where, formally, kn = 0). In such a state, the charge 
and spin densities are almost homogeneous along the z- 
direction. Furthermore, there are subgap states energet- 
ically located within the surface mode gap (A^). These 
states are localized near both caps and show interest- 
ing spin texture. Second, we observe significant out-of- 
surface components for the spin density associated with 
all energy eigenstates. Note that an out-of-plane spin tex- 
ture is on ly exp ected when trigonal warping effects are 
important ,1^2123 ggg also very recent experimental results 
reporting such features !^^'^^' However, to lowest order in 
momentum (around the F point), trigonal warping can be 
neglected, while we persistently find broken spin-surface 
locking also at the lowest energy scales. Our observations 
are instead related to the presence of non-differentiable 
sections of the surface, where the wavefunctions in the 
"trunk" and "cap" regions of the cylinder have to be 
matched. Such non-differentiable surface parts also ap- 
pear in the samples studied in Refs. 17181 and therefore 
surface probe experiments for these devices could di- 
rectly test our prediction of broken spin-surface locking. 
For diffcrentiable closed surfaces, we expect spin-surface 
locking to stay intact. In fact, the explicit solution of 
the pro blem in a spherical geometry exhibits spin-surface 
lo eking .fi^l^ We note in passing that for a flat TI surface, 
a time-dependent out-of-plane spin component can also 
be generated by elastic disorder. However, this compo- 
nent will precess around the momentum-dependent spin- 
orbit axis (which lies in the plane) and averages to zero on 
time scales corresponding to the inverse Fermi energy.l^ 

The structure of the remainder of this paper is as fol- 
lows. In Sec. ITTl we describe the results of our numerical 



calculations for cylindrical TI nanowire dots based on 
the effective low-energy theory. In Sec. |III[ an analytical 
approach based on the surface Dirac fcrmion picture is 
described for the same geometry. The spectrum and the 
spin density profile for the resulting eigenstates will be 
derived, and the results are compared to the numerical 
findings in Sec.|Tl] In Sec. |IV[ we turn to the tight-binding 
calculation and compare those computations to the pre- 
vious results. Finally, we conclude in Sec.|V] Appendix [A| 
contains a derivation of the surface Dirac fermion Hamil- 
tonian for an infinitely long nanowire using the approach 
of Ref. M 

II. EFFECTIVE LOW-ENERGY DESCRIPTION 

A. Model and numerical approach 

In this section, we compute the band structure of 
a cylindrical nanowire of length L and radius R fr om 
the effective low-energy theory of Zhang et oL'^SEII us- 
ing parameters for Bi2Se3. Up to terms of order with 
k = {kx, ky,k), the low-energy bulk Hamiltonian employ- 
ing the four bands energetically closest to the F point 



Hb = (kC^QTo + MfcO-QTz + [Ao{kxax + kyCjy 



C2k\ and Mk 



(2.1) 

Mo + Mifc2 + 



ky. This also defines the Fermi 



where €k — Cq + Cifc^ 
M2fc]^, with k\ = kl 

velocities vi — Bq/H (in z direction) and V2 = Aq/H (in 
the xy plane), see Eq. ( |l.ip . The basis states encoding 
the spin-parity structure, where Pauli matrices ct^ (r^) act 
in spin (parity) space and ctq (to) denotes the respective 
unit matrix, are explicitly given in Ref. 1211 This work 
also describes the extension of Eq. (2.1 1 to the case of 



eight bands or towards including trigonal warping. Using 
the parameters in Ref. [211 listed for convenience also in 
Table |lj the criterion for a TI phas^ is satisfied. One 
therefore must have an odd number of conducting surface 
modes when boundaries are present. 



Our aim is to describe the band structure of a finite- 
length cylindrical nanowire with axis along the z direc- 
tion. Due to rotational symmetry in the xy plane, it is 
useful to switch to cylindrical coordinates (r, 0, z). The 
"cylindrical" Pauli matri ces (g ,. , cr^ , ) then represent 
the physical spin operator,'2Sl211 



(2.2) 



and we refer to their local expectation values as "spin 
densities" below. The conserved total angular momen- 
tum operator is 

J = e~''^-'^/2 {-id^) e^'^-'^/2 ^ _^Q^ ^ ^j2. (2.3) 

For the finite-length cylinder we then construct the eigen- 
functions to the Hamiltonian (2.1) with Dirichlet bound- 
ary conditions, ^(r) = 0, on the surface, i.e., for \z\ < 
L/2 with r = R (cylinder trunk) and for \z\ = L/2 with 
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r < R (caps) . This is automatically achieved by expand- 
ing states in a complete orthonormal basis, {'ipai'r, 4>, z)}, 
that satisfies these boundary conditions. The quantum 
numbers a = (j, v, n, a) include the half-integer angu- 
lar momentum j, a radial index S N, the longitudinal 
quantum number 71 g N, and the spin index a — ±. Ex- 
plicitly, see also Ref. JJj for r < R and \z\ < L/2, the 
basis is chosen in the form 

MrA,z) = ^/^sin[vrn(z/£-l/2)]e^-^ -^7(^7"/^\ 

(2.4) 

where m = j — cr/2, V = irR^L is the cylinder volume, 
and jmu denot es th e vth zero of the Bessel function J^,. 



The basis set (2.4) satisfies the orthonormality relation 



JytPr '4'l{i')'4'a/Jr) = ^aa'- In addition, the basis states 
acquire a spinor structure in parity space not shown ex- 
plicitly in Eq. (lOl). 



Expanding the Hamiltonian Hi, [Eq. (2.1 1] in this basis, 
we obtain a matrix representation that allows for numeri- 
cal calculations in a truncated basis set. Upon increasing 
the basis set, numerical results for the spectrum turn out 
to converge rather slowly. We have performed a lattice 
regularization as in Ref. H7| in order to obtain manage- 
able matrix dimensions. Typically, we achieve conver- 
gence with « 8000 basis states for given j. The solution 
of the eigenvalue problem then yields the discrete en- 
ergy spectrum of such a quantum dot, E = Ej,s,±, where 
s G N labels the different states for the conduction or va- 
lence (±) band with given angular momentum j. Taking 
averages with respect to the corresponding eigenvector 
\'^j,s.±) then yields the spatially dependent charge den- 
sity profile for this state, {p){r,z). In addition, one ob- 
tains the local spin densities, {aa){r, z) with a = r, 0, z, 
and the local parity densities, {Tp)(r,z) with /3 — x,y,z. 
Rotational symmetry implies that all these averages are 
independent of the angular variable (j). 



B. Numerical results 

We now present the results of our numerical calcula- 
tions. The discrete energy spectrum for a TI nanowire 
dot with i? = 20 nm and L = 44 nm is shown in Fig. [T] 
The Kramers degeneracy results in an identical spectrum 
for j — >■ —j but with reversed spin and parity (r^) di- 
rections. We therefore show only the j > solutions in 
Fig.[T] Moreover, we focus on the topologically protected 
surface fermion modes inside the bulk gap A^. 

There are several noteworthy points about Fig. [T] 
First, comparison with our analytical results, see 



Eq. (3.191 and Sec. Ill below, shows that most levels are 



approximately recovered from the bulk dispersion rela- 
tion [Eq. (1.1 )] by simply imposing the standard quanti- 
zation condition /c„ = mr / L with n G N on the longitudi- 
nal momentum k. However, here additional states corre- 
sponding to n = emerge. These zero-momentum states 
are absent for Schrodinger fermions in a box. Note that 
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FIG. 1: (Color online) Energy spectrum Ej^s,± from numer- 
ical diagonalization of Eq. ( |2.1| for a cylindrical TI quantum 
dot with _R = 20 nm and L = 44 nm. Open blue triangles 
correspond to the numerical results, while filled red circles 
show the analytical prediction in Eq. ( |3.19[ ). Dashed lines in- 
dicate the surface gap for the respective angular momentum 
(j) mode in an infinitely long nanowire. The spin texture for 
the zero-momentum state with j — 1/2 indicated by the lower 
left arrow is shown in Fig. |2] The spin texture for the next 
higher state (higher left arrow) is shown in Fig. [s] The exam- 
ple of a subgap state (for j = 3/2, right arrowjis addressed 
in Fig.g 



for each j, there is precisely one n = state for the con- 
duction band and one for the valence band. Inspection of 
the density profiles for these states reveals almost homo- 
geneous charge, spin, and parity densities as a function 
of the z coordinate. Second, for j > 1/2, we find a pair 
of nearly degenerate subgap states inside the surface gap 
Ag. (The near-degeneracy is not visible in Fig. [l]) Such 
a subgap state is localized with equal occupation proba- 
bility at both caps. Furthermore, electron-hole symmetry 
is broken under the Zhang model Hh, in contrast to the 
analytical model in Sec. |III[ This is the main reason for 



the existing discrepancies between Eq. (3.191 and the nu- 
merical results, see Fig. [l] In fact, we have also carried 
out additional numerical calculations for an electron-hole 
symmetric version of Eq. ( |2.ip , where the corresponding 
results fit almost perfectly to the analytical results in 



Sec. |III| In particular, all subgap states then disappear. 

Inspection of the spin densities, {aa){r,z), and parity 
densities, {Tj3){r, z), for a given eigenstate {j, s, ±) yields 



(cr0)(r,z) = (T^)(r,z) = 0, 



(2.5) 



i.e., spin is never oriented in the circumferential direc- 
tion. In addition, there is now a finite radial spin com- 
ponent within the trunk region {\z\ < L/2) and a finite 
z-component within the caps (|z| = L/2). Hence, in gen- 
eral, the spin direction for a surface state points out of 
the surface: spin-surface locking is broken in this geome- 
try. This finding is in striking contrast to what happens 
in an infinite nanowire^^ and for a sphere.'^ Specifically, 
in the infinite cylinder case, the results corresponding to 
Fig. [2] show ((Tr) = refiecting spin-surface locking. 
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FIG. 2: (Color online) Spin density components (cr^) [(a)] 
and {(Tr) [(b)] in the rz plane, for the zero-momentum state 
with j — 1/2 indicated by the lower left arrow in Fig. [l] 



FIG. 3: (Color online) Same as Fig. [2| b ut for the next higher 
energy level (upper left arrow in Fig. nj) . 



Figure [2] shows the spin density profile for the lowest- 
lying ("zero-momentum") conduction band state with 
j = 1/2, indicated by the lower left arrow in Fig. [T] 
We indeed find an almost homogeneous spin density pro- 
file along the trunk, where spin is mostly aligned along 
the (negative) z-direction, see Fig. [2ja) . Also the charge 
density is practically homogeneous along the z-direction 
(data not shown). However, there is also a finite ra- 
dial spin component breaking spin-surface locking, see 
Fig. ^h). For the cap region, spin is mostly aligned 
along the radial direction, but again an out-of-plane com- 
ponent, now oriented along the z-axis, is clearly visible. 
For comparison. Fig. [3] shows the respective results for 
the next higher energy level (upper left arrow in Fig. [T]). 
Again we observe that spin-surface locking is violated, 
while in the infinite cylinder case one finds (cTj, = 0) 
(spin-surface locking) . 

For j > 1/2, our numerical results include an almost 
degenerate pair of subgap states, where the degeneracy 
is on top of the Kramers degeneracy. The charge density 
is then localized with equal probability near each of the 
two cylinder caps. A typical example for the spin texture 
of such a subgap state is shown (for j = 3/2) in Fig. |4] 
The out-of-plane spin part is identical on both caps, see 
Fig. 4 ^a) , but the in-plane (radial) component shown in 
Fig. 4 h) has opposite direction. 



Comparing our numerical results for parity, charge and 
spin densities, we find that for each eigenstate, they are 
linked by a set of general relations. In particular, for the 
(radially integrated) densities in the trunk region, we find 



{Ty){z) cx -{ar){z), {t^){z) cx (p)(z). 



(2.6) 



Note that in the infinite wire case^^^ the parity struc- 
ture is trivial in the sense that (for large R) the only 
non-zero component is (r^), cf. Appendix [A| 



III. SURFACE DIRAC FERMION THEORY 



In this section, we analyze the finite-length nanowire 
geometry of Sec.|TT]within a surface Dirac fermion theory, 
where we retain parity in the Hilbert space of the surface 
Hamiltonian. While this is not necessary for an infinite 
cylindei'^ or for the flat 2D surface,'^ the discussion at 
the end of Sec. |II B| shows that this extension is important 
here. 
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the results of Ref. [201 Evidently, both time-reversal and 
inversion symmetry are only kept intact when choosing 
T e {to, T^} in Eq. (3.1 1. We here set T = t^, as follows 



from the analytical derivation of Eq. (3.11 in Appendix 
1X1 as well as from numerical calculations based on the 
Zhang model for the infinite nanowire case.'^ 



Using Eq. (3.2), we can now switch to a ID represen- 
tation for a given angular momentum {j) channel. For 
given energy the ID spinor obeys the ID Dirac 
equation H^^^^pj 



Eipj wit 



]V2 



(3.4) 




1.0 



FIG. 4: (Color online) Same as Fig. [2] but for the "subgap" 
state with j = 3/2 indicated by the right arrow in Fig. [l] 



A. Infinitely long wire and symmetries 

We start with the parity-extended surface Dirac 
fermion Hamiltonian for an infinitely long nanowirej ^^ ^ ^^ l 



Hn = 



K 



r, 



(3.1) 



where the total angular momentum operator J has been 
defined in Eq. (2.3). T — acts in parity space and is 



determined below. We note that H]j respects all sym- 



metries present in Hh [Eq. (2.1)]. Specifically, azimuthal 
symmetry implies [HdjJ] — 0, and states are classified 
by half- integer j. 



V'((/),2:) = e 



e^^Ujiz), (3.2) 



jez+i/2 



with the ID spinor tpji^)- Time-reversal symmetry im- 
plies [Hjj,(d] = 0, with the time-reversal operator 8 = 
iffyToC, where C denotes complex conjugation. Finally, 
Eq. ( |3.1| exhibits inversion symmetry, [Hd,X] = 0, with 
the inversion operatoi^^SI 



Here, inverts the z coordinate, z - 



(3.3) 

-z, and Ra, shifts 



)-|-7r. The parity structure in Eq. (3.3) follows from 



see also the derivation leading to Eq. (A3) in Appendix 
\K\ Note that the representation of the Dirac matrices 7^ 
in Eq. ( A3 ) in terms of products of spin and parity ma- 



trices, see Eq. (3.4), is multi- valued and leads to a double 
counting of all surface states derived from Eq. ( 3.4 ) . Nev- 
ertheless, it is technically convenient to proceed in this 
representation, since the double counting can be easily 
circumvented, see Sec. [Ill CI 

The general solution to the ID Dirac equation with 
Hamiltonian H'^^^ in Eq. (3.4) reads 



V'j(z) 



-iA 



jikz 



-ikz 



, (3.5) 



with arbitrary complex coefficients {Ai, ^2 , -Bi , -82 ) • The 
spinors above are in parity space while Xj &cts in spin 
space. 



^. = cos 7, 
> i sin 7j 



3V2 

2RE 



3V2 

2RE 



(3.6) 



In Eq. (3.5), the longitudinal momentum k — k{E) fol- 
lows from vik — E"^ — {jv2/RY- Below we consider 
energies where k is real and positive. For a description 
of the subgap states discussed in Sec. IIIBI evanescent 
modes need to be studied instead. In the ID represen- 
tation, the inversion operator [Eq. (3.3)] becomes X 
with 



(3.7) 



Since Eq. (3.4) stays invariant under inversion. 



[ij(''\l] = 0, the eigenf unctions (3.5) can be classified 



as inversion symmetric or antisymmetric (cr = ±), 



Xi,^"\z)=a^,'P{^z). 



(3.8) 



From Eqs. (3.5) and (3.8), after a short calculation, we 



can therefore infer relations between the coefficients for 
given inversion symmetry (a): 



Ai = crBi 



A. 
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A general inversion-symmetric (antisymmetric) state 
thus takes the form 



(+) 



Ai cos7j cos(fcz) 

—iA2ax I I — sin7j sin(fc2:) 

iAi \ ( cos7i sin(fcz) 

I \ sm7j cos(A;zj 



(3.9) 



Both are parametrized by two complex numbers (Ai and 



A2), where the first (second) spinor refers to parity (spin) 
space. 



Next let us take the local expectation value for the 
charge, spin, and parity operators in a given general 
eigenstate i'^^'^ ■ First, with sin^(7_j) = cos^(7j), the 
charge density is 



(p)(z) = + |v42|') [cos2(7,,)cos2(fcz) + sin2(7,,)sin2(fcz)] . 

For the spin density, we obtain 

I -i(|Ai|2 + |A2|2)sin(27,)sin(2fcz) 

(^0) \{A = A 

(^3> / V (l^iP - 1^2!') [cos2(7,,)cos2(fcz) -sin2(7,,)sin^(fcz)] 

Similarly, the parity density is obtained in the form 

(t^) \ / crIm(AiA5)sin(27j)sin(2fez) 

\Ty) \{z)=\ crRe(AiA^)sin(27j)sin(2fcz) 
(r.) j \ - |A2|2) [cos2(7,,)cos2(fcz) +sin2(7,,)sin2(fcz)] 



(3.10) 



(3.11) 



(3.12) 



At this stage, the above results hold for an arbitrary 
inversion-symmetric (cr = +) or antisymmetric (cr = — ) 
state. Remarkably, the circumferentially oriented spin 
density {a^ always vanishes. This is in accordance with 
our numerical observations, see Eq. (2.5). The current 



density along the z direction can be obtained from the 
local operatoi!^^ 



(3.13) 



and therefore vanishes identically for these states as 
well.!^ This result stays valid for arbitrary inversion- 
symmetric boundary conditions (which do not mix V'j^'' 
states) . 

In order to reach agreement with the numerical results 



in Sec. II B the coefficients yli_2 should obey the three 
relations 



Im(AiA;)=0, Re(AiA;)>0, ^ IA2I. (3.14) 



B. Matching trunk and cap states 



The coefficients Ax 2 in Eq. (3.9 1 as well as the energy 
spectrum and the corresponding eigenstates can now be 
obtained analytically by matching the trunk states, see 
Eq. (3.91 for \z\ < L/2, with cap states at z = zLL/2. 



Each cap is described by a surface Dirac Hamiltonian of 
the form 



cap 



V2 



J 

(Tr H 



(3.15) 



The parity matrix is uniquely determined by imposing 
time-reversal and inversion symmetry, and also appears 
in Eq. (2.1). In the angular momentum (ID) represen- 
tation, for given energy E, a cap state ^'^'^''(r, z — C,L/2) 
[with C = i and inversion symmetry index a — ±] takes 
the general form (for a — C, — —, this should be multi- 
plied by —1) 

^Y\r, Ci/2) ^C^[\^® 0(r) + C_c {\) ® 0(0 

(3.16) 



Indeed, the first relation implies consistency with with complex coefficients C± and the spin spinor 
Eq. (2.5). The second relation ensures that {Ty){z) oc 



{ar){z), see Eq. (2.6). The third relation is required to 
have non- vanishing spin-parity densities (cr^) and {tz){z). 
Moreover, notice that then {tz){z) cx {p){z), in accor- 



J,-x/2mrlv2) 
isg-a{E)Jj+il2{\E\r/v2) 



(3.17) 



dance with Eq. (2.6) 



The other spinors in Eq. (3.16) are in parity space 
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Matching the trunk states [Eq. ( 3.9 1] and the cap states 
[Eq. ( |3.16 l] for given j and a by continuity at z = ±L/2 
and r = R,we obtain four Hnear equations for the four co- 
efhcients (^i, A2, C+, C_). A nontrivial solution follows 
when the corresponding determinant vanishes, yielding 
the condition {u = \E\R/v2) 

(jUM^) + J'+iM^)) sin(27,) sm{kL) = 0. (3.18) 

For real-valued fc > 0, this equation can only be satisfied 
when sin(/cL) = 0. This implies the standard longitu- 
dinal momentum quantization condition kn{E) — mr/L 
with n S Nq. The corresponding eigenenergies then fol- 
low from the bulk dispersion relation in Eq. (1.1), 



are also reproduced, see Eq. (2.6). Specifically, the dot 



K- 



±^{nnvi/LY + {3V2/Rf. (3.19) 



For a given level, the wave function amplitudes A12 and 
C± then satisfy three conditions plus the overall normal- 
ization constraint. With p — o'(— )" = ± and the above 
definition of m, we have (7+ ~ pC- and the relation 

J.j+p/2{u)Ai +p sgn(i;) Jj_p/2iu)A2 = 0. (3.20) 

Moreover, for p = + the third condition reads 

cos(7,)Ai - 2(-i)"Jj_i/2(u)C+ = 0, (3.21) 

while for p = — this instead becomes 

sin(7,)Ai - 2(-^)" sgn(S) Jj+i/2(u)C+ = 0. 

These relations determine all possible wavefunctions for 
the closed cylinder surface. 

Of course, since we considered only solutions of 



Eq. (3.18) with real fc, subgap states were not captured. 



However, within the linear-in-fc approximation underly- 
ing the approach here, we find that there are no subgap 
states at all, i.e., the corresponding matching problem 
with evanescent trunk modes does not permit a nontrivial 
solution. The numerical approach of Sec.|TT]also indicates 
that in order to obtain subgap states, it is necessary to 
include higher-order terms (in fc) breaking electron-hole 
symmetry in the Hamiltonian. 



C. Spectrum and eigenstates of the dot 

Using the above explicit solution for the wavefunction 
of the complete cylinder, we obtain the coordinate de- 
pendence of all densities of interest. We here describe 
the results for the trunk region only. First of all, we 
recover 



see Eq. (2.5). Moreover, the relations 

(Ty) (X -{C7r), (r^) OC (p) 



eigenstates have the energy E = Ej^n,± specified in 



Eq. (3.19). Wc show below that these energies are not 



degenerate, i.e., only one specific inversion parity a given 
by 



o- = aj.n^± = T(-l)"sgn(j) 



(3.22) 



will be physically realized. With fc = nn/L > 0, we find 
for the charge and spin densities {u — \E\R/v2) 



{p)(z) cx 1 - sgn(i;)— cos(2fcz), 
u 



{(Tr){z) CX -v/l — p /v?' sin(2fcz), 



(3.23) 



{az){z) cx cos(2fcz) 



-sgn{E). 



These results are in good agreement with the numerical 
results for the spin texture obtained for the Zhang model 
in Sec. |ll] In particular, they show explicitly that spin- 
surface locking is broken. 

Interestingly, for each total angular momentum j, 
there are two zero-momentum states corresponding to 
conduction and valence band, respectively. Their inver- 
sion symmetry properties are determined by 



cr = -sgn(j£;j^„=o,±), 



(3.24) 



since for states with the opposite value of a, all densities 
in Eq. (3.23) vanish. For the physically allowed fc = 
state with a in Eq. ( |3.24 ), from Eq. (3.23) we instead find 
spatially uniform densities {p){z), (r^) (z), and {az){z), 
while all remaining spin or parity density components 
vanish. 



We now compare the densities in Eq. (3.23) to the 
numerical results in Sec. [llj and also address the dou- 



ble counting problem mentioned in Sec. Ill A Eq. (3.20) 



shows that indeed Im(yli742) = 0, in accordance with our 
numerical results in Sec. IIB see Eq. (3.14). Also the 
relation \Ai\ 7^ IA2I in Eq. (3.14) is evidently satisfied. 



However, not all the states can be realized physically, as 
is clear by comparing to the condition Re(AiA2) > in 
Eq. (3.14) found numerically in Sec.|llj In order to under- 
stand this restriction, we note that the operator S = cr^r^ 
commutes both with the cap Hamiltonian [Eq. ( |3.15 )] 
and with the inversion operator I [Eq. (3.7)]. This im- 
plies by continuity that trunk states [Eq. (3.9)] at the 
end points {z = C^/2 with ( — zL) are eigenstates of 5 
as well, 



2^S(Ci/2) 



p^(;>(C£/2), 



where the eigenvalues p = (t(— 1)" follow from Eq. (3.9 
and the definition of S. For n — 0, however, Eq. (3.24 
implies that only the eigenvalue p — — sgn(j_Ej^o,±) is 
physically realized. By continuity, this value must also 
apply for the full Hilbert space of conduction or valence 
surface bands. We therefore obtain the "selection rule" in 



Eq. (3.22) restricting the Hilbert space of allowed states. 
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This explains the condition Re(AiA2) > ^ind resolves 
the double-counting problem. We mention in passing 
that the latter problem is automatically avoided when 
retaining terms of order fc^ in the Hamiltonian, where 
the spin-parity eigenstates with p — ± have different en- 
ergy. In the Dirac theory, this implies a "spontaneously 
broken symmetry" encoded by Eq. (3.22 1. 



D. Effective boundary conditions 



It is also possible to derive the results in Sec. |III C| 
without explicit construction of the cap states. To that 
end, let us briefly consider a class of general boundary 
conditions at the cylinder ends, z = with C = ±, 

by imposing the local gauge constraints 



V^(0,Ci/2)-AcV('/',Ci/2), 



(3.25) 



where A^; = A^^ = A^. The spin-parity structure of A^ 
can be determined by requiring time-reversal invariance, 
[A^,0] = 0, and invariance under inversion, [A,^,!] = 0. 
In addition, we require the boundary operator to com- 
mute with Hd, which is a natural assumption for closed 
surfaces. As a result, with arbitrary angles ri±, we find 

A^ sin(?y^)a-oTz + cos(?7^)a-rTy. 

Passing to the ID representation, i.e., for given half- 
integer j, these constraints read 

= AcV',(CL/2), (3.26) 
Af = sm{ri^)aoT^ +cos{7]^)axTy. 

Applying the inversi on op erator I, see Eq. (|3.7[), to the 



boundary condition ( |3.26[ ), we find A_|_ = A_ and hence 
r]± = rj. Only then will inversion symmetry be pre- 
served for the confined states. The parameter r] (with 
< 77 < tt) cannot be fixed by symmetry considera- 
tions alone but depends on the physical boundary con- 
dition imposed at the ends, i.e., the boundary matrix A 
effectively encodes the matching of trunk states with cap 
states. Contrary to the commonly employed boundary 
conditions ,'^233 ^j^g operator A commutes with the cur- 
rent operator jz, see Eq. (3.131, while the anticommuta- 



tor is always nonzero. Since the boundary conditions are 
invariant with respect to inversion, they do not mix the 
states (3.9) with opposite inversion parity a. Using the 



boundary condition (3.26), some algebra yields for both 



solutions and for arbitrary energy E the condition 



cos r/ 
1 — sin 7] 



-A2 = 0. 



(3.27) 



Comparing this to Eq. (3.20), the energy-dependent an- 



gle 77 can be explicitly related to the above wavefunc- 
tion matching procedure, and the subsequent results in 
Sec. |III C] can be obtained under a purely ID description 
of the trunk states alone. 



IV. MICROSCOPIC TIGHT-BINDING 
APPROACH 

A simple microscopic model for a strongTI was pre- 
viously proposed by Fu, Kane, and Mele.'^S' The model 
consists of a single-band tight-binding model on a di- 
amond lattice and includes spin-orbit couplings. With 
lattice fermion operators Ci (spin is kept implicit), this 
Hamiltonian has the form 



[4 



x4]) 



where a is the cubic lattice constant, tij are hopping pa- 
rameters connecting nearest neighbors, and the last term 
describes spin-orbit coupling of strength Ago through a 
second-neighbor hopping between sites which de- 
pends on the two nearest-neighbor vectors d}'^ connect- 
ing those two sites. In order to generate a full gap in the 
bulk spectrum, a distortion tj^ — » t + 5t\s introduced for 



dij along the (111) direction.'*^ 

To define the nanowire, we proceed as in Ref. [121 by se- 
lecting the growth direction ez along the (111) axis and 
keeping all sites within a given radius R. The unit cell 
of the infinite nanowire thus defined contains six planes 
of sites corresponding to the three stackings of the two 
fee sublattices of the diamond lattice, and has the pe- 
riod dcoii = •s/Sa. Although this model is not completely 
equivalent to the one of Eq. (2.1), quantitatively similar 



behavior of the surface gap as a function of R has been 
obtainecP^ by setting a = 2.8 nm and —2t + 5t = Mq. 
The finite dot geometry is then defined by setting the 
length of the nanowire along the (111) direction to a 
given value L. To maintain the aspect ratio of the cylin- 
drical dot studied numerically in Sec. |llj we set R = 3a 
and L = 4-\/3a. This corresponds to a cluster of 1592 
sites, which is a small enough size to keep a reasonable 
computational cost of the calculations while still allowing 
for a meaningful comparison of the spin texture with the 
results of Sees. HIl and IIIII 

The band structure of the infinite nanowire and the en- 
ergy levels in the finite dot geometry are shown in Fig. [5] 
where we again focus on states energetically inside the 
bulk gap. As expected, both the bands of the infinite wire 
and all dot levels are twofold Kramers degenerate. For 
the chosen parameter set, we find two subgap states ap- 
pearing inside the surface gap A^. (Since there is no full 
rotational symmetry anymore, we cannot classify states 
by j here.) 

To get an idea of the spin texture of the dot states, we 
now focus on the two states indicated by arrows in Fig. [5] 
The lower arrow corresponds to a subgap state, and the 
second one refers to the lowest-lying state within the con- 
duction band of the infinite wire. The corresponding spin 
densities are shown in Figs. |6] and [Tj As in the effective 
low-energy theory and the surface Dirac fermion descrip- 
tion, we again observe (cr^) = 0, and therefore only (cTz) 
and (ar) are shown. Note that in the infinite wire case, 
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FIG. 5; (Color online) Band structure of the infinite TI 
nanowire (red solid curves) and energy levels (denoted by 
open black triangles) of a finit e dot . These results follow from 
the tight-binding model [Eq. ( |4.1[ )] with wire axis along the 
(111) direction. The dot length is set to L = 4%/3ffl, and the 
radius is J? = 3a. The arrows indicate the two levels whose 
spin texture is analyzed in Figs. |6] and [T] respectively. 



one would instead find (ct^) = consistent with spin- 
surface locking.^^ 

The subgap states have a charge density mostly lo- 
calized near the caps of the cylindrical dot, again with 
out-of-plane (in-plane) spin components that are identi- 
cal (oppositely directed) on both sides, reproducing the 
results of Sees. |ll] and |III| However, unlike the continu- 
ous model, the tight-binding model predicts a spin tex- 
ture with a superimposed atomic-scale oscillation. This 
oscillation stems from the finite k — Tr/dcoii value at the 
Dirac point in this model. 

On the other hand, for the lowest-lying state within 
the conduction band, corresponding to the "zero- 
momentum" state in Sees. and |III[ Figure [7] shows that 
the density is largest along the cylinder trunk. The spin 
is predominantly oriented along the negative z-direction, 
but with a finite oscillatory component in the radial di- 
rection that breaks spin-surface locking. 

To conclude, even though the results of the tight- 
binding model are not fully equivalent to the ones in 
Sees. |TT] and |III[ we find that the main properties of 
the subgap and lowest conduction band states are re- 
produced. 



V. CONCLUSIONS 

In this paper we have studied the band structure of 
a quantum dot made of a strong topological insulator 
using three different approaches, namely the low-energy 
theory of Zhang et q/, i^ofzi] g^^^ effective surface Dirac 





i.a 



FIG. 6: (Color online) Spin texture in the rz plane obtained 
from the tight-binding model [Eq. ( |4.1[ )]. The spin densities 
{(Tz) [(a)] and (ar) [(b)] are shown for the subgap state indi- 
cated by the lower arrow in Fig. [S] 



fermion theory, and numerical calculations for a tight- 
binding model on a diamond lattice with strong spin- 
orbit couplings. The considered geometry, with flat caps 
terminating a finite-length cylindrical nanowire, is char- 
acterized by sharp edges where the cylinder trunk and 
caps meet. Such edges are also present in typical "meso- 
scopic" TI devices studied experimentally.'^KI^ All three 
approaches show that spin-surface locking is generally 
violated due to presence of these edges. As also found 
in a recent ab initio study,'^ a finite reflection proba- 
bility for Dirac fermions in each part results when two 
surfaces are patched together. In our case, we have a 
Fabry-Perot-like setup where standing waves can build 
up. The resulting spin density then exhibits spatial os- 
cillations reminiscent of a spin density wave state. The 
spin direction of the oscillatory parts points out of the 
surface while non-oscillatory spin density contributions 
stay locked to the surface. 

The spectrum of such a quantum dot shows sev- 
eral surprising features. First, when starting from the 
band structure of the infinitely long wire [Eq. (1.1)], 



imposing the usual longitudinal quantization condition. 
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1.0 




FIG. 7: (Color online) Same as Fig. |6] but for the state indi- 
cated by the higher arrow in Fig. [5] 



kn = nn/L, here allows for a nontrivial eigenstate with 
fc = n = 0. In such a zero-momentum state the charge 
and spin densities along the trunk are basically homo- 
geneous. Second, albeit the surface bands of an infinite 
nanowire exhibit a gap, the finite-length nanowire dot 
has subgap states when electron-hole symmetry is bro- 
ken. The wavefunction of such a subgap state is local- 
ized on both caps simultaneously. The obtained energy 
spectrum and corresponding spin textures are important 
ingredients for a theory of mesoscopic transport through 
TI dots. In general, we also expect Coulomb interactions 
to be relevant, in particular charging effects should be 
visible. We plan to address these questions in the future. 
Moreover, extensions of the theory to include an applied 
magnetic field, where the typically large and anisotropic 
Lande factop^H implies that the Zeeman field is crucial, 
are also left for future work. 
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TABLE II: Representation of spin-parity matrices in terms 
of Pauli matrices 7^ acting in surface-state subspace. 



Appendix A: Surface Hamiltonian for infinite 
nanowire 



Here we analytically derive Eq. (3.1) for an infinite 
nanowire using the gap-inversion model of Ref. 1181 Our 
starting point is Eq. (2.1) to linear order in k (we put 
Co = here), 



Hb = McToTz + [vika^ + V2{kx(Jx + kyay) 



(Al) 



Following Ref. UHl we assume that the gap parameter 
changes sign at r = i?, i.e., M{r) = Mq sgn(i? — r) with 
Mq > 0. For r < R {r > R), the material is then in 
the topologically nontrivial (trivial) phase. In the ID 
representation, for given j and k, solutions to the Dirac 



equation take the form e e 



E,e'^''^V,(0, with 



the ID radial Dirac equation 
Here k = corresponds to 



A/(r)o-oTz + V2 



1 

2^ 



For fc = 0, there are two solutions 



IS cos 7 
— sgn(A/) sin 7 



where the first (second) spinor refers to parity (spin) 

space, sin 7 = \J\- cos 7 = \J\ + T^-i ^^^d V2K. = 

\J Mq — E'^ . For r > R, we have to replace Ij-1/2 ^ 
Kj-1/2 and Ij+1/2 — >■ —Kj+i/2 in the spin part, where 
I,^ and Ki, are modified Bessel functions. The general 
solution is then given by 



V'j(r) 



aiV]+^(r)+/3iV^j-)(r), r < R, 
a24'l'^Hr) + M^l'\r), r > R, 



(A2) 



where the coefficients ai 2 and /3i_2 are obtained by re- 
quiring continuity of the wave function at r = i?. This 
results in a linear system of equations for the coefficients, 
which has a nontrivial solution under the condition 

{cos'^{l)Ij-i/2K.j+i/2 - sin^(7)/,-+i/2iv:j_i/2) 
X (sin2(7)/j_ 1/2^4:^+1/2 - cos^ {-/)!. 



J+1/2KJ-1/2) = 0, 
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where all Bessel functions have the argument kR and 7 
is evaluated for M = Mq. Assuming kR ^ 1 and us- 
ing the asymptotic form of the Bcssel functions, we find 
cos(27) ~ ±j/{K.R) and thus £'j,±(fc = 0) = ±jv2/R. 
The corresponding k = Q wavefunctions, il>j,±{f'), are 
given by Eq. (|A2|) with j3i — ±ai, j32 — ±Q!2, and 



ai 



-kR 



0-2 =T 



^kR 



27r 



For small A: 7^ 0, the effective surface Hamiltonian, 



Eq. (3.11, is obtained by projecting Eq. (Al) onto the 



subspace spanned by the above fc = states. We thereby 
find 



JV2 



(A3) 



where the 7^ are Pauli matrices in the zero-momentum 
subspace. In this way, all combinations of spin-parity 
matrices, aiTj, can be represented in the truncated basis. 
For kR ^ 1, we obtain the results in Table [III Of course, 
this representation is not single- valued, i.e., there is no 
one-to-one correspondence between the aiTj and the 7^ 
matrices. In particular, the spin operators ay^^ a-nd the 
parity operator are entangled, since CFy^zTo — and 
CToTz = 0. Taking % = OyT^ and 7^ = -ct^t^, cf. the 
last column of Table |ll) one gets th e sur face Hamilto- 
nian quoted in the main text, see Eq. (3.4). Because this 



representation is multi-valued, the replacement of 7^^^ 
by products of spin and parity matrices causes a double- 
counting of states if used naively. This double-counting 
problem is, however, not severe and can be circumvented, 
as we discuss in Sec. IlllCl in the main text. 
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